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We investigate the surface critical behavior of two-dimensional multilayered aperiodic Ising models 
in the extreme anisotropic limit. The system under consideration is obtained by piling up two types 
of layers with respectively p and q spin rows coupled via nearest neighbor interactions Ar and A, 
where the succession of layers follows an aperiodic sequence generated through substitution rules. 
Far away from the critical regime, the correlation length ^x, measured in lattice spacing units, in 
the direction perpendicular to the layers is smaller than the first layer width and the system exhibits 
the usual behavior of an ordinary surface transition. In the other hmit, in the neighborhood of the 
critical point, diverges and the fluctuations are sensitive to the non-periodic structure of the 
system so that the critical behavior is governed by a new fixed point. We determine the critical 
exponent associated to the surface magnetization at the aperiodic critical point and show that the 
expected crossover between the two regimes is well described by a scaling function. From numerical 
calculations, the parallel correlation length is then found to behave with an anisotropy exponent 
z which depends on the aperiodic modulation and the layer widths. 
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I. INTRODUCTION 

Since the discovery of quasi-crystalsJ3 the critical prop- 
erties of aperiodic systems have been intensively studied 
(for a review, see Ref. ||). These systems display, in some 
cases, the same type of singularities which can be encoun- 
tered in random systems, the study of which is a quite 
active field of research. 

In layered systems, generalizing the McCoy-Wu 
model,Hi3 aperiodic distributions of the exchange inter- 
actions between auccessive layers in the Ising model have 
been considered. Oila In Ref. O, the Onsager logarithmic 
singularity of the specific heat was found toJpe washed 
out. The major result was obtained by Luckjlj by a gen- 
eralization to layerecL^erturbations of the Harris criterion 
for random systems.tJ 

According to Luck's criterion, aperiodic modulations 
may be relevant, marginal or irrelevant, depending on the 
correlation length exponent u of the unperturbed system, 
and on a wandering exponent uo which characterizes-the. 
fluctuations of the couplings around their average.tatll 
Systematic studies of the surface critical properties for ir- 
relevant, marginal, and relevant perturbaticps have been 
achieved in the extreme anisotropic limit .EScil In the case 
of bulk marginal sequences (i.e. which modifies both sur- 
face and bulk properties), an anisotropic scaling behavior 
was obtained in the Hamiltonian limit ,E2rEi as well as in 
the classical twQ-dimensional counterpart ,cfl or in hierar- 
chical models.cJ In these systems, both surface and bulk 
properties have been studied and the bulk specific heat 
exponent, given by a = 1 — z, where z is the anisotropy 
exponent, is always negative leading to a non-diverging 
singular behavior. 

In the present paper, we continue the investigation 
of the surface critical behavior of marginal sequences, 



through a study of the crossover between usual ordi- 
nary surface transition and aperiodic critical behavior 
in multilayered Ising models. Up to now, only layered 
systems have been considered, and we study here mul- 
tilayered systems made of the piling up of two kinds of 
layers of widths p and q. In such systems, together with 
the perpendicular correlation length the first layer 
width p plays the role of a relevant length scale entering 
the description of the critical behavior in the vicinity of 
the free surface. Depending on the value of the scaled 
variable p/^±, one thus expects a crossover between the 
two regimes. In section O, we recall some generalities 
about aperiodic sequences generated through substitu- 
tion and give a short summary of the relevance criterion 
proposed by Luck. The critical behavior of the surface 
magnetization in multilayered m odel s at the aperiodic 
fixed point is determined in Sec. and the crossover 
effect, which constitutes the main result of the paper is 
presented in Sec. IV. In section 0, a numerical study 
of the parallel correlation length ^|| allows the determi- 
nation of the anisotropy exponent of multilayered Ising 
models in agreement with a relation, already established 
in the layered case, which seemingly has a general field 
of application. 



II. LUCK'S CRITERION 

Many results have already been obtained in 2d layered 
Ising models with constant interactions Ki along the lay- 
ers and modulated interactions K2{k) between successive 
spin rows, k and k+\. The transfer matrix should in prin- 
ciple be diagonalizcd for arbitrary couplings ,E3 but in the 
following one considers the extreme anisotropic limit. In 
the Hamiltonian limit, Ki oo, K2 0, while keeping 
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the ratio = K2 (k) / K* fixed, where = — ^ In tanh Ki 
is the dual coupHng, the row-to-row transfer operator 
T = exp[— 2if J'TY] involves the Hpniltonian of a quantum 
Ising chain in a transverse field :lIIl3 



n 



1 ^ 

^EH + ^^-^fe'^^V) ' (2.1) 



k=l 



where the cr's are Pauli spin operators. j—. 

The Jordan- Wigner transformation pj f ollowed by a 
canonical transfoKcuation, maps Eq. ( |2.l[ ) onto a free 
fermions problemn3 



(2.2) 



The fermion excitations £c 
the linear system 



follow from the solution of 



eatf'a(fc) = -<?c(fc) - Afe<?c.(fc + 1) 
£a^a{k) = -\k-l^a{k - 1) - ^a{k) 



(2.3) 



which may be rewritten as a single eigenvalue equation, 
either for <I> or 4" which are assumed to be normalized 
and free boundary conditions Aq = Al = are imposed. 

The surface properties are defined as usually by ma- 
trix elements: the surface magnetization is given 
by (l|crf |0) where |0) is the ground state of TL and |1) 
is the first excited state. Using the transformation to 
?7-fermionic operators, the surface magnetization takes 
the simple form nis = <Pi{l). The surface energy den- 
sity Cs involves the lowest two-fermion excited state 
2) ~ ViVtl^)^ ^iid can thus be written = (2|crJ|0) = 
(£2 -ei)^'i(l)<?2(l). 

With an aperiodic modulation of the couplings, one 
may write: 



At 



(2.4) 



where fk, which may be or 1, is determined by the 
aperiodic sequence. Let — J2k=i fk be the number 
of modified couplings on a chain with length L; their 
asymptotic density is poo = limi^oo '^■l/-^ The critical 
coupling Ac obeys the following conditionEJ 



lim n(AOc 



/L 
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which, using equation (2.4), gives: 



A. 



(2.5) 



(2.6) 



The sequences considered below are generated through 
substitutions rules, either on digits ("period-doubling") 
or on pairs of digits ("paper-folding"). ____ 

The "period-doubling" (PD) sequenccBBEl'Elmay be 
defined by the substitutions 5(0) = 11 and S{1) = 10, 
which, starting on 1, give successively: 



n = 1 1 

n = 2 10 

n=3 10 11 

n=4 10 1110 10 



(2.7) 



The "paper-folding" (PF) sequence is generated 
through substitutions on pairs of digits, and details can 
be found in Ref. |23| , |33| . 

Most of the propertiestjoLa sequence are obtained from 



its substitution matrixllj'EZi with entries giving the 
number nf^''^ of digits (or pairs) of type i in S{j). The 
asymptotic density of modified interactions, poo is de- 
duced from the eigenvector corresponding to the leading 
eigenvalue Ai of the substitution matrix, al low ing the 
calculation of the critical coupling via Eq. (2^). The 
length of the sequence after n iterations of the substitu- 
tion rules is asymptotically related to the leading eigen- 
value through Ln A". Finally, the fluctuations of the 
/fe's at a length scale Ln are governed by the next-to- 
leading eigenvalue A2: 
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(2.8) 



where F{x) is a periodic "fractal function"0 as shown in 
Fig. 0. 
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FIG. 1. Fluctuations of the sum of the fk's around 
their average for three different aperiodic sequences (top: 
Rudin-Shapiro, uj = 1/2; middle: "paper-folding", uj = 0; 
bottom: Fibonacci, ui = —1). The symbols correspond to 
sizes which do respect the symmetries of the chains and the 
dotted lines have slopes uj. 

The cumulated deviation per bond of the couplings 
from their average A is thus written 



1 L„ 

'5(in)--|-E(^fc-^)^'^|A2| 



5L-^ 



(2.9) 



fe=i 



where 6 = A(r — 1) is the amplitude of the modulation 
and UJ is the wandering exponent given by: 
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(2.10) 



Since the relevant length scale is the correlation length 
^, the aperiodicity thus induces a thermal perturbation 



t 



1 + v{lj-1) 



(2.11) 



involving a crossover exponent (f)^ When (f> > 0, the per- 
turbation diverges as t —>■ 0, corresponding to a relevant 
perturbation. When </) < 0, the perturbation is washed 
out when the critical point is approached; it is irrelevant. 
When (/) = 0, the perturbation is marginal and may lead to 
a non-universal behaviour with perturbation-dependent 
critical exponents. In the following, we will only consider 
marginal sequences. 

In the case of the 2d Ising model, since v ^ 1, the 
relevance of the perturbation is directly determined by lu 
as shown in Fig. [|. 



III. SURFACE MAGNETIZATION IN 
MULTILAYERED ISING MODELS 

A. Critical point and surface magnetization 

The first excited state of the Hamiltonian, |1}, is de- 
generate with the ground state in the ordered phase of 
the i nfin ite system and, since a simplification occurs in 
Eq. ( |2.3| ) when ei vanishes, the first relation provides a 
recursion for the components of $i. After normalization 
one obtain^he surface magnetization of the semi-infinite 
system as:Eil 



5-1/2^ S- 



j=i fc=i 



(3.1) 



which remains valid for any distribution of the interac- 
tions. 

We consider here a multilayered Ising model where the 
aperiodic sequence fk — 1, 0, 1, 1... refers now to suc- 
cessive layers with p, q, p, p . . . spin rows. A layer la- 
belled by index k contains p (respectively q) exchange 
interactions Xr (respectively A) between nearest neigh- 
bor rows according to the following scheme showing the 
first three layers: 



p terms 



Xr Xr 



.Xr 



k 

fk 
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p terms 



Xr Xr 



...Xr 

.2p + q 
(3.2) 



Let us first determine the critical coupling Ac. The 
width of layer k is written 



mk = {p- q)fk + q- 



(3.3) 



On the system containing L layers, one has layers of 
type p and L — of typ e q and the critical value is still 
determined through Eq. ( p.5| ) where 



l[(^i)'J'' - (Ac 
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In the thermodynamic limit, it yields 



Ac = r (p-9)poo+<? 



(3.4) 



(3.5) 



The aperiodic series (3T) entering the surface magne- 
tization has the form 
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^M(A,r) = l + ^(Ar-^0"" + (Ar-^0"""'I](Ar-^^) 
1=1 1=1 

-H(Ar/i)-2™i(Ar/^)-2™2 ^(Ar^')-2' + . . . (3.6) 



1=1 



where the superscript (m) is for the multilayered system. 
This expression may be shortened with the notations 



o-fc = ^(A? 



(Ar)2-1 



1=1 



ftp if fk = 1 

if fk - 

(3.7) 



and gk = (Ar^*" ) 2™'' . Equation ( ^^ ) then becomes 

(A, r) = 1 + 6,E(A, r) + (a^ - 6,)E'(A, r) (3.. 

where 

00 j 

S(A,r) = ^n5fe, 

j=0 k=0 



E'(A,r) = ^/,+in5fc, 

j=0 k=0 



(3.9) 



and /o = 0. In the second sum, S', /j+i may be rewritten 
as 



fj+i - 



(AP-9rP)-2 _ 1 
and Eq. (|3.8|) translates into 



(3.10) 



<jMm ,A - 1 _ A2'?(ap - bq) 



- l)(ap - b,) 



(XP-irP)-^ - 1 



E(A,r). (3.11) 



The critical behavior of the surface magnetization is 
governed by the sum E, for which an identity can be 
found, relating it to the aperiodic series S'-^^ of the usual 
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layer ed s ystem in the case p = q = 1. Taking account of 
Eq. (2^), the series S'(') (A, r) in Eq. takes the form 



2p 



(3.12) 



where the upperscript (I) stands now for the layered prob- 
lem. Ort&, uS | yia lly studies this expression for specific se- 
quences Jl^'EllE^ but for our purpose a more convenient 
expression can be rewrittei|U-as a function of the devia- 
tion from the critical point :E3 



t = l- 



Equation (3.12) becomes 



5(')(i,Ac) = ^(l-i)-'A-20"-«./p~). 



(3.13) 



(3.14) 



The analysis of the sum S is now facilitated, since the 
following identity holds: 



E(A,r) = S^^'>{X'',XP-''rP). 
With the variable t, it becomes 



(3.15) 



S(i, Ac) = ^(1 - t)(P-'?)"^+«(A«)^2(^'~"^/''~' (3.16) 

and the critical point behavior reduces to the same series 
as the layered problem up to the transformation Ac ^ A^: 

oo 

S(0, Ac) = ^(A^)-^(^'-"^/''°=) = S^^\0, XI). (3.17) 

This expression is the basis of further asymptotic. analy- 
ses using the scaling method proposed by Igloi.Lj Let a 
power series M{x) have a power law singularity at the 
transition point x = xq- M{x) ~ (1 — x/xq)^ . The trun- 
cated series of the first L terms at the critical point thus 
have the following behavior Ml{xq) ^ L^^ which allows 
the determination of [3. Equation ( 3.17] ) also holds for 
the truncated series and leads to a simple relation be- 
tween the critical exponents of the layered problem and 
the multilayered one: 



/3(™)(Ac)=/3y^(A?) 



(3.18) 



B. Period-doubling and paper- folding sequences 

The period-doubling sequence, given after several sub- 
stitutions in Eq. ( ^.7| ), corresponds to a vanishing wan- 
dering exponent and an asymptotic density of modified 
couplings poo = 2/3, leading to the critical co uplin g for 
the multilayered Ising model couphng via Eq. (3.5): 



(3.19) 



Equation (|3.16|) can be rewritten as an infinite product ,E£l 
whose first / terms contain, at the criti cal point t — Q, 
the first L — 2^' terms of the sum ( 3.14 ) which becomes 



(0, Ac) = [(1 + Xf){\ + A-^)]' ^ (22')2/3<'"' . (3.20) 
The surface exponent then follows :!! 

Going back to the original parameters, one finally ob- 
tains: 



/3^)(r,p,g) 



1 / ^'^ \ 

infr^P + Q -\- T 2p-|-q j 

21n2 ■ 



(3.22) 



The surface magnetization exponent depends on the am- 
plitude of the coupling ratio r through Ac as expected for 
a marginal sequence, but also on the layer widths. The 
variation of /3s™'' for several values of the parameters is 
shown in Fig. 0. 




FIG. 2. Marginal variation of the surface magnetization ex- 
ponent /3s"''' with the coupling ratio r for the period-doubling 
multilayer. The singularity is always weaker than in the ho- 
mogeneous system and the surface transition is of second or- 
der for any value of r,p, q. 



The analysis of the series ( p. 16 ) cas, be completed 
by a usual finite-size-scaling argument ,E3 which implies 
that the following critical point behavior also holds: 
Sl(0, Ac) ^ L^^"". The anomalous scaling dimension 
of the surface magnetization Xm^ = Ps/v then takes the 
same value than j3s, and this requirement imposes that 
the correlation length exponent in the perturbed system 
keeps its unperturbed value v = 1. This property holds 
for any marginal sequence and can be understood within 
Luck's criterion, which has to remain valid in the per- 
turbed system: the wandering exponent being constant, 
4> — Q indeed implies the constancy of v. 
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In the case of the paper-folding sequence with po 
1/2, the critical point is determined by 



Xr = r P+1 . 



(3.23) 
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FIG. 3. Marginal variation of the surface magnetization 
exponent Z?!'"^ with the coupling ratio r for the paper-folding 
multilayer. The singularity is weakened by the perturbation 
when r < 1. 

The same scaling analysis has been done in the case 
of the paper-folding sequence for the layered problem.c^ 
Translated to the multilayered problem, it leads to 



Sl=2<(0,A,) = (1 + A2?)' 
and to the surface exponent: 



(2 



2 In 2 



In 1 



2pq_ 

p+q 



2 In 2 



(3.24) 



(3.25) 



whose variations are shown in Fig. ^. While period- 
doubling is a symmetric sequence where the relation 



Ps'{r,p,q) = (3s'{r~^,p,q) holds, it is not the case 
of the paper-folding perturbation for which the surface 
magnetic exponent has a monotonous variation. The sin- 
gularity is thus weakened for r < 1 only. This result can 
be understood if one looks to the density of perturbed 
layers (of type p). Their density close to the surface is 
indeed larger than the asymptotic density Poo , leading to 
a weaker average coupling in the vicinity of the surface 
than in the bulk when r < 1. In the other regime r > 1, 
the average coupling close to the surface is enhanced and 
a stronger singularity follows. 



IV. CROSSOVER IN MULTILAYERED ISING 
MODELS 

The study of aperiodic multilayers does not present 
more difficulties than the usual layered systems, but it 



can produce an interesting crossover in the critical be- 
havior. While the perpendicular correlation length is 
of the order of the size of the chain in the close neighbor- 
hood of the critical point, in which case the critical prop- 
erties are governed by the aperiodic fixed point, this is no 
longer true when the system is moved far away from Ac. 
In this latter situation, the correlation length decreases 
as A increases and, at some point, becomes of the order 
of the first layer width p. The aperiodic structure of the 
system then becomes irrelevant and the behavior is con- 
trolled by the semi- infinite homogeneous fixed point, the 
unperturbed exponent /3s pi 1/2 of the ordinary surface 
transition being recovered.EJ 

The evidence of this phenomenon appears immediatly 
on a log-log plot as shown in Fig. ^ in the case of period- 
doubling sequence. The values of r and p — q have been 
chosen for clearness in order to keep a constant value 
= 2. 

The surface magnetization can thus be written as fol- 
lows: 



m,{t,p)^e^^f{p/U) 



(4.1) 



where f{x) is a scaling function involving the relevant 
length scales which describe the surface behavior. It is 
shown on a log-log scale in Fig. ||where m^t"^/^ is plotted 
against pt and ~ t~^. 

For a fixed value of /Si™', one obtains two different 
universal curves, one for r > 1 and the other for r < 
1. This is due to the properties of the period-doubling 
sequence which allows to reach the same M"^ > 1/2 
with two different values of the coupling ratios r at fixed 
p and q values. The paper-folding sequence exhibits a 
rather different behavior shown in Fig. 0. 




FIG. 4. Crossover between homogeneous fixed point 
{f3s = 1/2) towards aperiodic fixed point (here /Js"*' = 2) 
when the critical point is approached (period-doubling se- 
quence). The layer sizes are p = g = 10 (r = 2.295), p — 20 
(r = 1.515), p = 50 (r = 1.181), p = 100 (r = 1.087), p = 200 
(r = 1.042), p = 500 (r = 1.017), p = 1000 (r = 1.008), 
p = 2000 {r = 1.004) and p = 5000 (r = 1.002). 
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FIG. 5. Universal curves for the surface magnetization 
in the period-doubhng sequence with Z?!™' = 2 and 1 and 
p = (? = 10, 20, 50, 100, 200, 500, 1000, 2000, 5000. 



The scaling function of Eq. (4.1) has to satisfy appro- 
priate asymptotic behaviors in order to provide a correct 
description of the two limiting power-law behaviors of 
771s. Far away from f = it has to keep a constant value 
(the usual surface magnetization amplitude of the homo- 
geneous system) to agree with 777^ ^ t^/^. This require- 
ment is satisfied in Fig. || and |^ and is responsible for 
the vanishing slopes for large values of pt. On the other 
hand, in order to recover the right varying exponent /Ss™'* 
close to t = 0, the following power-law is expected in the 
vicinity of the critical point 



(m) 



1/2 



(4.2) 



The average slopes calculated numerically in this second 
regime are in agreement with the expression of 9 given 
above. 




FIG. 6. Universal curves for the surface magnetization 



in the paper-folding sequence with f3. 



(m) 



3,2,0.4,0.2 and 



q = 10, 20, 50, 100, 200, 500, 1000, 2000, 5000. 



The scaling functions shown in Figs. || and || exhibit a 
quite good agreement with the required asymptotic be- 
haviors. In the homogeneous surface regime (far away 
from t — 0), this fact should be surprising since t is 
the deviation from the aperiodic critical point, given in 
Eq. ( 3.13 ). Introducing the homogeneous critical cou- 



pling Ac.o = 7" of the first (perturbed) layer, and 



1 — (Ar) ^, one has 



^A^ is a measure of 



the corresponding deviation to 
t ~ to{l + 5/to), where S = 1 ~ r 
the difference between the two critical point locations. 
Since 5 is a small quantity in our calculations (for exam- 
ple S ~ 0.005 in the case p = q = 1000, r = 1.008 in 
Fig. 1^), we can expand Eq. (4.1) up to linear order in S: 



ms{t,p) ^ t, 



1/2 



1 



2to) ^ {u 



(4.3) 



and there appears simply a correction to scaling term 
which is proportional to the difference between the two 
critical couplings. With the range of values used here, 
this correction to scaling term is negligible. 

We can also point out log-periodic oscillations exhib-. 
ited by the scaling function in the "aperiodic" regime.ESI 
This behavior is related to the discrete scale invariance 
of aperiodic sequences. The oscillating amplitude. A, 
is usually a periodic function of lnS^j_/ Inm where 777 is 
the discrete scale factor (4 for period-doubling and 2 for 
paper- folding). In this multilayered system, the correla- 
tion length is naturally measured in units of p such that 
the scaling function in the neighborhood of the critical 
point can be rewritten 



f{x) ^ A(hix/ \iim)x^ 



leading to 



In /(a;) ^ ln[yl(lna;/ln 777)] +61tix. 



(4.4) 



(4.5) 



in agreement with the oscillations around an average 
slope 9. 



V. CORRELATION LENGTH AND MARGINAL 
ANISOTROPY 

In Ref. the low-energy excitation behavior has been 
studied for the layered problem, leading to a finite-size 
power law at the critical point involving an anisotropy 
exponent z. The conjecture 

was proposed,^ and it has then been proved for specifie. 
sequences, using renormalization group calculations .cll 
and has furthermore been supported for an arbitrary 
marginal perturbation by an approximate calculation of 
ei(L)Jl3 leading to ei{L) ^ ■ms{L)rhs{L) YlfS^ X^^ from 
which Eq. (5.1) follows. Moreover, the constancy of the 



6 



correlation exponent in the perturbed systems has been 
obtained exactly in Ref. 

Assuming the same finite-size behavior for the first 
component of the eigenvectors corresponding to the two 
lowest excitations, $i(l) ~ $2(1) ^ the surface 

energy density exponent also follows 



(5.2) 



The correspondence between the classical two- 
dimensional system and the quantum chain shows that 
the correlation length in the Euclidean time direction 
(along the layers) is given by the inverse of the small- 
est non- vanishing gap in the Hamiltonian spectrum, i.e. 
^11 = l/si at Ac- In the low-temperature phase A > Ac, 
the lowest excitation ei vanishes exponentially due to the 
asymptotic degeneracy of the ferromagnetic ground state 
in the thermodynamic limit, and the parallel correlation 
length is thus given by the next fermion excitation: 



£2 



(5.3) 



In the previous sections, we found the marginal behav- 
ior of the surface magnetization at the aperiodic critical 
point, and its crossover towards the usual ordinary sur- 
face transition fixed point in the ordered phase. We can 
expect here, as in the usual layered problem, a strong 
anisotropic scaling behaviour responsible for a special be- 
havior of at the aperiodic fixed point. It can be studied 
through the temperature dependence when A > Ac {via 
£2) or by finite-size scaling at the critical point A Ac 
(via £1). 



With free boundary conditions, the system (2.3) can 
be rewritten in a matrix form {A — B)(A + -B^q, = 
e^$Q. Here, the standard notation has been used.EJ Due 
to the tridiagonal structure of the "excitation matrix" 
{A — B){A + B), the recutaion relations for the compo- 
nents of (f>cv take the formO 



<?«(/ + 1) 



„ - 1) 







1 



Si ti{£a) 



(5.4) 



with 

A/-1 , , 

Sl = r — , ti(ea) 



and boundary terms 



2<l <L-1 
(5.5) 



>a(l) 

XML + 1) 



Ti 



<?a(0) 



(5.6) 



<?a(0) = 

Ai<?>c.(L + l) = 0. 



Finally, one has 



<^a(i) 
/(£) 



■T2T1 



(5.7) 



(5.8) 



where the eigenvectors are not normalized at this step 
and the excitation energies are deduced from the positive 
zeroes of /(e). In a multilayered structure, within a given 
layer of size p, the sam e m atrix appears p ^ I times and 
in the case p — q, Eq. (5.8) becomes 



■■■Qfu.iPhr-'Qf.jAPhr-'Ti (5.9) 

where P/^ is a "propagation" matrix inside layer k and 
^ fk+i fk ^ "transfer matrix" from layer k to layer k + 
1. We can easily give a closed form for the propagation 
terms: 



1 

-1 tk{£) 



after diagonalization, leads to 



fk 











Jk 



(5.10) 



(5.11) 



where a± are the eigenvalues of P and S the chang- 
ing of basis matrix. 




IpLe/n 

FIG. 7. Typical shape of the function /(e) for pe- 
riod-doubling sequence, with L = 2'^,p = q = 2, r — 1.1, 
and several values of A in the ordered phase. 

It is now possible to get Ml for any value of the layer 
sizes p, and, in the case of large values of p, this technique 
is quite efficient compared to a direct diagonalization of 
the excitation matrix. Some numerical troubles can nev- 
ertheless occur, since the function /(e) exhibits oscilla- 
tions whose amplitude increases sharply with the sizes L 
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and p while the period decreases rapidly. An exemple is 
shown in Fig. for a small size. It is clear from Fig. 
that, as A increases from its critical value, the lowest exci- 
tation El decreases while the next one, £2, increases. The 
temperatm'e behaviour of i^n is in agreement with a mod- 
ified power-law ^|| ~ t~^, where the anisotropy exponent 
z is assumed to be the sum 

z{r, p, q) = (r, p, q) + {r-\ p, q) (5.12) 

as in the layered system. The temperature dependence 
of e^^ is shown on a log-log scale in Fig. || in the ordered 
phase A > Ac. From the linear behavior, we can obtain 
the anisotropy exponent z. The numerical results de- 
duced from the slopes of the asymptotic linear behavior 
are given i n Ta b. ^. The numerical data are in agreement 
with Eq. (^.12|) . 

s I 1 1 1 1 1 




2 - 



-4.5 -4 -3.5 -3 -2.5 -2 -1.5 

\nt 

FIG. 8. Temperature dependence of the parallel correlation 
length on a log- log scale in the ordered phase A > Ac. The 
deviation from the asymptotic linear behavior (dashed lines) 
when t — » is a finite size effect. 



TABLE I. Parallel correlation length exponent //y , deduced 
from the slopes of log-log plots of v.s. t, compared to the 
anisotropy exponent z for the period-doubling sequence. 



r 


p = q 




z 


1.1 


10 


1.102 


1.072 


1.1 


20 


1.286 


1.273 


1.1 


30 


1.493 


1.575 


0.9 


10 


1.099 


1.087 


0.9 


20 


1.334 


1.330 


0.9 


30 


1.626 


1.686 



values of L (2^° for the period doubling sequence for ex- 
ample) and extrapolation to infinite size has been per- 
formed, using the BST algotrithm, which is efficient for 
finite-size-scaling analyses.CJ 



TABLE II. Anisotropy, surface magnetization and surface 
energy density exponents deduced from a finite-size-scaling 
study for the period-doubling sequence. 



r 


p = q 


2 

num. 


theor. 


X 

num. 


m) 

theor. 


(m 

num. 


theor. 


2 


2 


1.15(2) 


1.149 


.58(3) 


.574 


2.29(1) 


2.298 


2 


3 


1.33(2) 


1.322 


.69(5) 


.661 


2.62(3) 


2.644 


3 


2 


1.36(2) 


1.357 


.70(4) 


.678 


2.70(2) 


2.713 


3 


3 


1.75(2) 


1.737 


.90(4) 


.869 


3.45(2) 


3.474 



VI. CONCLUSION 

We have presented the study of the influence of layer 
widths p and q in aperiodic multilayered Ising mod- 
els. The perturbations under consideration correspond 
to vanishing wandering exponents, and thus lead to 
marginal critical behaviors in the two-dimensional Ising 
model. From a finite-size scaling analysis of the ape- 
riodic series which defines the surface magnetization, we 
obtained the corresponding critical exponents which con- 
tinuously vary with the modulation ratio. This defines 
the aperiodic fixed point behavior characterized by an 
anisotropic scaling which has been studied numerically. 

When the system, in the ordered phase, is moved away 
from this critical point, the surface magnetization ex- 
hibits a crossover towards another regime, governed by 
the ordinary surface transition fixed point. We have 
shown that the crossover is well described by scaling func- 
tions of a single scaled variable which obey well-defined 
asymptotic behaviors. 

Aperiodic multilayers are now feasible by molecular 
beam epitaxy,Lj and although nothing experimental has 
been done up to now to study the critical behavior of 
such systems, crossover effects between the two regimes 
should be experimentally observable. 

One should finally mention that the same type of study 
for relevant aperiodic modulations would be interesting 
especially when an enhancement of the modified cou- 
plings close to the surface leads to first order surface 
transitions in a two-dimensional system. This will be 
the subject of a further investigation. 



A finite-size-scaling study provides an alternate numer- 
ical calculation of these exponents. 

It has been done in Tab. || for the anisotropy expo- 
nent, but also for the surface magnetization and en ergy 
density exponents, in order to check the conjecture (5.2) 
for the multilayered system. The layer widths have been 
limited to small values in order to reach large enough 
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